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ABSTRACT 

Evolution of fast, radiative shocks in high density medium is presented. Ionizing spec- 
tra and approximate broad band light curves of the shocked gas are calculated. Emer- 
gent shock spectra, as seen by a distant observer, are obtained from photoionization 
models. The emergent spectra have a power-law shape F v oc v~ a with mean spectral 
index a ~ 0.6 — 1.0 in the energy range 0.01 — 10 keV, and have a high-energy cut- 
off corresponding to the original shock velocity. It is shown that the models exhibit 
promising features that may account for some photometric and spectral properties of 
Active Galactic Nuclei. 

Key words: hydrodynamics - shock waves - instabilities - ISM: supernova remnants 
- galaxies: active - X-rays: galaxies 



1 INTRODUCTION 

In the series of papers Terlevich and his group (see, e.g., 
Terlevich et al. 1992) developed the starburst model for Ac- 
tive Galactic Nuclei (AGNs). In this model AGN is powered 
by compact, dense supernova remnants (cSNRs), and the 
bulk of radiation is emitted by the supernova shock wave 
evolvin g in a dense medium (no ~ 10 7 cm -3 ). Terlevich 
et al. (1992) calculated ID and 2D hydrodynamical mod- 



els of cSNR evolution and demonstrated that it was possi- 
ble to recover observed characteristics of Broad Line Region 
of AGNs. Using very simple assumptions Terlevich et al. 
(1994a, TTRFM) successfully explained observed difference 
between moments of maximum continuum and line emis- 
sion taking into account dependence on ionization parame- 
ter. While the rapid X-ray variability of AGNs still remains 
t o be explained in this model, preliminary investigations 
(Terlevich et al. 1994b) showed that interactions between 



dense, fast moving clumps of gas could be partially responsi- 
ble for observed variations of high-energy emission. Detailed 
critici sm of the starbur st mo del was presented by Heckman 
( |l99l[ ) and Filippenko ( |l992j ). 

Apart from their hypothetical presence in AGNs, com- 
pact supernova remnants were detected in nearby galaxies. 
On the basis of spectroscopi c studies suc h objects like SN 
1980K ( |Fesen fc Becker 1990[ ), SN 1987F ( |Filippenko 1989| ), 
or SN 1988Z (Stathakis & Sadler 1991; Turatto et al. 1993) 
are interpreted as supernova remnants evolving in a medium 
of density as high as 10 7 cm -3 . The lack of appropriate the- 



oretical models was noted by Leibundgut et al. (1991) in 



their study of SN 1986J and SN 1980K. Both the possible 
role of cSNRs in AGNs and the remark by Leibundgut et al. 
gave motivation to our work on hydrodynamical evolution 
of fast shocks evolving in a dense environment. 

The sub ject i s not new. Over a decade ago Chevalier 
& Imamura (1982) showed that radiative, steady shocks are 
subject to an oscillatory instability. This result was con- 
firmed on the basis of nonlinear hydrodynamical analysis 
for nonstationary shocks (Gaetz, Edgar & Chevalier 1988; 
Stone & Norman 1993) as well as for steady radiative shocks 
(Innes, Giddings & Falle 1987). Both groups found unstable 
behavior of shocks faster than ~130 km s . The spectra of 
faster shocks (v s < 1100 km s -1 ) were discussed in some de- 
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tail by Binette, Dopita & Tuohy ( |1985| ). They found that the 
emergent shock spectrum can be approximated by a power- 
law with a spectral index a ~ 0.5 for energies lower than 
~ 1 keV. However, all theoretical investigations mentioned 
above apply to relatively slow shocks in low density medium 
(n < 100 cm -3 ). Therefore, they are not applicable to con- 
ditions expected for cSNRs or for central regions of AGNs. 

In the first paper of the series (Plewa & Rozyczka 1992, 
Paper I) we investigated the evolution of fast (v s q > 1000 
km s -1 ) radiative shocks in a uniform medium of density 
n — 10 7 cm -3 . The shock and its post-shock region were 
found to be highly nonstationary, leading to the conclusion 
that some stochastic component might be present in the 
evolutionary pattern, but we failed to obtain any definite 
results due to the low resolution of the models. Subsequently, 
we used a high-resolution adaptive grid technique (Plewa 
1993, Paper II) and we found quasiperiodic patterns in the 
evolution of the shocks. The flow behind the main shock 
was found to be highly discontinuous due to the presence 
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of secondary shocks leading to a stochastic behavior of the 
shock luminosity. 

In Paper I and II our model shocks were formed by 
the reflection of a steady stream of gas from a rigid wall, 
and they were not directly applicable to any astrophysical 
situation. In the present paper we improve the models by 
allowing for a secular decrease in stream velocity to mimic 
the evolution of the forward SN shock. Also, approximate 
spectra of shocks are obtained under the assumption of col- 
lisional equilibrium. Description of the method is given in 
Section ^[ In Section ^| and ^ we present details of the shock 
evolution, the broad band photometry and the spectrum of 
the radiation emitted by the shock, as well as results ob- 
tained from photoionization models. Possible relation of the 
models to AGNs is discussed in Section 5l 



2 THE METHOD 

The hydrodynamical evolution of the shock is modeled using 
Piecew ise Parabolic Method (PPM) of Colella & Woodward 
(1984) in ID plane-parallel geometry. Our preliminary in- 
vestigations presented in Paper I showed that even the high 
resolution of the PPM scheme is not sufficient to model ac- 
curately the most luminous and dense regions of the shocked 
gas if a standard uniform grid is used. The evolution of 
such regions is characterized by cooling time scale much 
shorter than the dynamical scale of the flow. For example, 
to avoid complications arising from the lack of resolut ion in 
strongly radiating regions Blondin, Fryxell & Konigl (199C) 
used a 



treatment lowering artificially local emission 
rates. Our numerical experiments showed that such pro- 
cedure strongly affects the results, and we concluded that 
to model the evolution correctly a much higher resolution 
should be used. 

We found that to resolve thin, strongly radiating regions 
of shocks faster than ~ 1000 km s _1 the grid should con- 
tain about 10 4 equally spaced zones. The cheapest way to 
obtain such a high resolution is the adaptive grid method. 
In our approach presented in Paper II the grid motion is 
governed by an adaptive grid algorithm developed by Dorfi 
& Drury ( 1987 ). The grid points are distributed accordingly 
to a user-defined resolution function. For our problem the 
resolution function was defined as the sum of gradients of 
density, internal energy, and emission rate (see Paper II for 
details). We also added the term proportional to the gas 
density to obtain higher resolution in high density regions. 
The grid equation was solved implicitly at the beginning of 
every hydrodynamical time step. The grid was divided into 
500 zones and we used Courant number equal to 0.8. The 
number of time steps needed to complete the evolution using 
maximum point concentration equal to 2.2 x 10 4 was about 
2 x 10 s . 

At the left boundary of the grid we impose a reflecting 
boundary condition, while from the right side the grid is fed 
with gas of constant density and temperature. To mimic the 
decrease of shock velocity in a way appropriate for SN shock 
we vary the velocity of the inflowing gas with time as 

_ 5 

Win =uoX (t/t Bg + 1) ~ , 

where i sg has the interpretation of the time that has to elapse 



Figure 1. Schematic representation of initial conditions for shock 
evolution. 



since the SN explosion before cooling effects become notice- 
able. According to Shull (1980) 



i sg = 0.63£sn (^r) 



-3/4 



yr, 



if the explosion energy is _Esn x 10 5 erg and the density of 
the ambient medium is no. This way we account for the fact 
that in ~ t sg years after the SN explosion the mean velocity 
of the forward SN shock begins to decline like (t/t sg )~ 5 ^ 7 
(Shull 1980; see also Fig. 3 in Falle 1981). 

This set of initial and boundary conditions leads to the 
formation of a strong shock at the left boundary of the 
grid (Fig. ^]). The initial velocity of the shock for constant 
gamma-law gas with 7 = 5/3 is given by v s o = |wo- Let us 
note that, when compared to real SNR, the left boundary of 
the grid can be regarded as the interface between the ejecta 
and the gas shocked by the forward supernova shock. Thus, 
the advantage of our approach is the possibility to concen- 
trate all grid points in the most interesting narrow region of 
the remnant where strong cooling occurs. The disadvantage 
is the simplified boundary condition imposed at the left edge 
of the grid, which in general may influence the structure and 
dynamics of the dense shell. This, however, does not hap- 
pen in the case of shock with v B o > 3000 km s _1 which form 
their shells far away from the left edge of the grid (Fig. ^). 
The reflecting boundary condition affects only the models 
of shocks with v s o < 3000 km s _1 , which form two shells, 
one of them being attached to the left edge of the grid (Fig. 
|§]a). However, as the latter shell is much less massive than 
the shell formed close to the main shock, the influence of the 
boundary condition on both the dynamics and luminosity of 
the post-shock region as a whole is marginal. 

Following Terlevich et al. (1992) and TTRFM, we as- 
sume for the ambient medium a density no = 10 7 cm' 3 
and a temperature To = 10 4 K. These values are admit- 
tedly arbitrary. However, the choice of temperature is practi- 
cally unrestricted because of high shock velocities (the shock 
would remain very strong even if To was an order of mag- 
nitude higher). As for the density, it was shown in Papers 
I and II that its reduction or increase results respectively 
only in increase or reduction of the evolutionary time scale 
and the size of the post-shock region, while all characteristic 
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evolutionary features of the model remain unchanged. The 
adopted value of no = 10 7 is particularly interesting because 
of indications that it may be characteristic of the medium 
surrounding cSNRs. 

Radiative losses are allowed for in the post-shock re- 
gion, while the shock moves to the right through the yet un- 
shocked gas. The losses are calculated implicitly assuming 
collisional equilibrium conditions for optically thin plasma 
with solar abundances. We use cooling function obt ained 
with the help of CLOUDY 84.06 code ( [Ferland 199S| ). We 
assume that the heating/cooling balance is achieved at a 
temperature equal to the temperature of the unshocked am- 
bient medium T m i n — To — 10 4 K. Therefore, in the regions 
where temperature has dropped below T m i n (e.g. due to gas 
expansion) the temperature was increased to reach the lower 
limit given above. 

Also under the assumption of collisional equilibrium, 
spectra of the radiation emitted by the hot shocked gas are 
calculated at selected evolutionary moments. These spectra 
will be called ionizing spectra in order to distinguish them 
from the emergent shock spectra obtained from photoioniza- 
tion calculations presented in Section 4.2 below (the latter 
are the spectra emerging from our models which would be 
registered by a distant observer if the space between the 
observer and the shock would be empty.) To save the CPU 
time, in place of the extremely time consuming CLOUDY 
code we use the following approximations for calculation of 
the ionizing spectrum. Continuum emission consists of the 
f-f emissio n wit h approximate Gaunt factor calculated after 
Hummer (1988), while the b-f and 2-photon emissions are 



calcul ated u sing formulae given by Mewe, Lemen & van den 
Oord ( 1986) ) . Contribution from over 2600 lines is calculated 
using line emissivit ies ta ken after Stern, Wang & Bowyer 
( [l978| ) and Mewe ( 199§ ) (see also Mewe, Gronenschild & 
van den Oord 1985). 

In order to estimate the relative broad band distribu- 
tion of energy in the ionizing spectrum we define a set of 
three filters corresponding to wavelength intervals A < 100 
A (high energy, HE), 100 < A < 905 A (medium energy, 
ME), and above 905 A (low energy, LE). The filter trans- 
mission function is defined as a fraction of the total energy 
emitted at a given temperature per given wavelength inter- 
val (LE, ME, HE; Fig. ^). The transmission functions are 
calculated by detailed integration of the spectrum emitted 
by the optically thin plasma as given by the CLOUDY code. 

Crucial to our approach is the assumption of collisional 
equilibrium which makes the problem numerically tractable. 
Admittedly this assumption does not work well in the lower 
temperature range. However, it is justified for the bulk of 
the shocked gas with temperatures higher than 5 x 10 5 K 
(only below that value the cooling function begins to be his- 
tory dependent; see, e.g., Shapiro & Moore 1976; Gaetz et 
al. 198 8). Moreover, as it was demonstrated by Gaetz et al. 
(1988) details of cooling are not very important for global 



shock dynamics (obviously, they play a more significant role 
in the evolution of cold regions with T < 10 5 K) . The cool re- 
gions are the source of ultraviolet and optical radiation only 
and we do not expect significant changes of the high energy 
emission when nonequilibrium effects would be considered. 
These issues are further discussed in Sect. [11]. We do not 
calculate preionization structure of the gas upstream from 
the shock. Also, we give here only results for solar composi- 



Figure 2. Definition of the broad band filters. HE, ME, and LE 
filter transmissions are drawn by thick, thin and dotted lines, 
respectively. 



tion, although higher metallicity would be more appropriate 
for both cSNR and AGN applications. Taking all the simpli- 
fications into account, our shock spectra are to be treated as 
first approximation only, especially in LE and ME ranges. 
On the other hand, we are fairly confident that the dynam- 
ics of the shocks has been modeled properly, and that any 
improvements in the treatment of cooling processes would 
not affect dynamical evolution of the shocks in a significant 
way. 



3 DYNAMICAL EVOLUTION 

We calculated models for initial shock velocities equal to 
1000, 2000, 4000, and 6000 km s" 1 . It should be stressed here 
that only the fastest shock model may be directly related to 
cSNRs evolving in a uniform medium of no = 10 7 cm -3 , 
as, according to Shull (1980), in such a case the velocity of 
the SN shock at t — t ss is equal to ~ 6000 km s _1 . The 
remaining models have no direct application, and they only 
serve to illustrate the strong dependency of the evolution on 
v s o (one may expect, however, slower SN shocks to appear 
in nonuniform media, e.g., in dense shells of wind bubbles 
blown by progenitors.) 

Here we present details of the evolution for shocks of 
velocities equal to 2000 (the slow shock) and 6000 km s _1 
(the fast shock), as the evolution in the two remaining cases 
shows no qualitatively new features. The evolution was fol- 
lowed up to evolutionary time t = 3 yr and t = 12 yr for 
slow and fast shock, respectively. Summary of the results for 
all the cases considered is presented in Table [l]. Note that 
the total shock luminosity was calculated assuming that the 
radius of the emitting sphere is equal to R s ^ = 3 x 10 16 cm 
(see Fig. 8a of Terlevich et al. 1992). Evolution of the model 
shocks in form of a sequence of vertically shifted density 
plots in logarithmic scale is shown in Fig. |^ (all plots were 
taken at equidistant intervals of time). The left boundary 
of the plot corresponds to the contact discontinuity between 
the ejecta and the ambient medium, and the ambient gas 
enters the grid from the right side. 
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Figure 3. Density evolution of model shocks, (a) initial shock velocity: v 3 q = 2000 km s — ; distance scale: 1.1 X 10 15 cm; time scale: 3 
years; density scale: 1.78 X 10 -18 — 2.90 X 10 -13 g cm -3 ; (b) initial shock velocity: v s o = 6000 km s — x ; distance scale: 1.2 X 10 16 cm; 
time scale: 12 years; density scale: 1.79 X 10 — 18 — 1.26 X 10 — 13 g cm -3 . The vertical scale may be read from the jump in the lowermost 
which corresponds to an increase of density by a factor of 4. See text for details. 



Table 1. Basic parameters of the model shocks. 



^s0 


^max 


^max 


^max 


km s — 1 


cm 


erg cm -2 s — 1 


erg s — 1 


1000 


1.8 X 10 14 


7.5 x 10 6 


9.4 x 10 40 


2000 


1.0 x 10 15 


3.7 x 10 7 


4.7 x 10 41 


4000 


4.7 x 10 15 


8.7 x 10 7 


1.1 x 10 42 


6000 


1.2 x 10 16 


1.5 x 10 8 


1.9 x 10 42 



ti s o - the initial shock velocity; iJ max - the maximum extension 
of the shock; F max — the maximum energy flux; L max — total 
shock luminosity assuming the radius of the emitting sphere 
equal to R s ^ = 3 X 10 16 cm. 



3.1 Evolutionary phases 

In all cases listed in Table |l| the evolution of the shock could 
be divided into three distinctive phases: 

(1) Nearly adiabatic shock expansion. The evolution starts with 
a nearly adiabatic shock expansion as cooling is ineffective 
at temperatures over 10 7 K for typical cooling function (see 
Raymond, Cox & Smith 1976). As both temperature of ra- 
diating gas and velocity of inflowing gas decrease with time, 
the shock slows down while cooling efficiency begins to grow 
with temperatures falling below 10 7 K. 

(2) Thin shell formation. At temperatures lower than ~ 10 6 
K the cooling time scale suddenly drops down. This phe- 
nomenon is known as catastrophic cooling (Falle 1975, 1981) 
and it results in formation of thin transition regions (here- 
after referred to as cooling waves) in which thermal energy of 
the gas is rapidly radiated away at nearly constant density. 
The cooling waves are the source of almost all the emission 
generated by the model in ME and LE energy bands. As 



our resolution function contains the term proportional to 
the luminosity gradient, cooling waves are well resolved. 

Soon after the formation of cooling waves the cool gas be- 
gins to condensate into a very thin and dense shell, visible as 
the highest peak in both panels of Fig. |^ (the less prominent 
spikes in Fig. ^ originate from minor flow discontinuities like 
weak shocks). At this moment the grid points strongly con- 
centrate in the region of shell formation (the maximum grid 
concentration corresponds to 2.2 x 10 4 equally spaced zones). 
This grid concentration is close to a practical upper limit, 
as at significantly higher resolutions the time step becomes 
prohibitively short due to the CFL condition. This phase 
of the evolution really challenges our adaptive grid method 
which has to redistribute the grid points between the main 
shock profile, the fast moving cooling waves, and the rapidly 
growing shell. All these structures are well resolved, proving 
the efficiency of the method. 

(3) Oscillatory instability of the main shock. Finally, loss of the 
pressure in the post-shock region due to cooling leads to 
the recession of the main shock, and to its reflection from 
the dense shell marked by a kink in shock position in Fig. 

in which position and velocity of the shocks are plotted 
against time. In Paper I we showed that the velocity of the 
reflected shock depends on momentary structure of the shell 
and parameters of the inflowing gas. At the moment of the 
shock reflection the velocity of the inflowing gas is by a factor 
of few lower than at the beginning of the evolution. This 
implies lower velocities of the reflected shocks (v IS < 1000 
km s" 1 ) and, therefore, shorter evolutionary time scales. The 
evolutionary time scale of the main shock is additionally 
reduced due to strong cooling in the region just behind the 
shock, as after the reflection the temperature of the post- 
shock gas never exceeds 10 7 K. The mean period of the 
oscillations is roughly equal to P osc w 14.5 d and P osc w 
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Figure 4. Temporal variations of position and velocity of the main shock (dashed and solid lines, respectively), (a) initial shock velocity 
Vso = 2000 km s _1 ; distance scale: 1.04 X 10 15 cm; velocity scale: 2000 km s _1 . (b) initial shock velocity d s o = 6000 km s _1 ; distance 
scale: 1.18 X 10 16 cm; velocity scale: 6000 km s _1 . 



3.3 d for fast and slow shock, respectively. The period of 
the main shock oscillations corresponds to the fundamental 
mode of the oscil latory instability discus sed by Chevalier 
& Imamura ( 1982 ) and Gaetz et al. ( 198E ). The oscillatory 
phase is preceded by a relatively short period of time during 
which the density contrast between the inflowing gas and the 
shell is too low for the shock reflection to occur (the reflected 
shock is very weak and its cooling time is extremely short). 
This phase is marked by a chaotic behavior of luminosity and 
velocity of the main shock. Semi-regular oscillations begin 
only after the density contrast has become high enough for 
an efficient shock reflection. 



3.2 Dependence of the evolution on shock velocity 

In both cases described above the evolution starts with a 
nearly adiabatic expansion of the main shock followed by 
formation of the thin shell. In the case of the fast shock 
(v s0 = 6000 km s" 1 , Figs. |(b), §(b), and |J>) the shell con- 
denses far away from the rigid left boundary of the grid, and 
it moves outwards due to high pressure in the hot, rarefied 
gas left behind it. This region (called hot cavity by Terlevich 
et al. 1992) is created as the post-shock temperature drops 
faster due to rapidly decreasing inflow velocity than due to 
cooling (as a result, catastrophic cooling occurs closer to the 
shock). This is not true for shocks slower than ~ 3000 km 
s" 1 as the initial post-shock temperatures are lower by an 
order of magnitude and, therefore, the cooling is much more 
effective. In that case the shell recedes to the left, and we 
observe the formation of the second shell at the left bound- 
ary of the grid. The collision of two shells at the end of the 
run for slower shock (upper left part of Fig. ^a) is marked by 
prominent luminosity spike at time t « 2.5 yr (Fig. ^,). Fi- 
nally even for shocks of higher velocities the shell will collide 
with the left boundary of the grid. 



PHOTOMETRIC AND SPECTRAL 
FEATURES 



4.1 Broad band emission 

The behavior of the broad band emission of ionizing radia- 
tion produced in the post-shock region is shown in Fig. ^[ 
During the first evolutionary phase associated with nearly 
adiabatic shock expansion the total luminosity smoothly 
increases, and the shock slows down. At the moment of 
shell formation LE/ME luminosities rapidly increase while 
HE emission remains nearly constant. Fast variability of 
LE/ME emission during shell formation indicates nonuni- 
formity of the gas entering the shell. This nonuniformity of 
the cool, condensing gas arises from small, numerical fluc- 
tuations present in the post-shock region of nearly standing 
main shock (see Fig. ^ which are amplified by catastrophic 
cooling. 

After shell formation (at t as 1 and t m 3.2 yr for slow 
and fast shock, respectively) the HE luminosity steadily de- 
clines. The decline rate is higher for the slow shock, as in 
this case the hot region behind the shell has lower tempera- 
ture. In the case of fast shock this hot, rarefied gas radiates 
the rest of its thermal energy at t » 10.5 yr. 

The main shock reflection from the shell (t ~ 1.6 yr, 
Figs, ^(a) and ||(a); t tti 5.0 yr, Figs. [I](b) and||b) is marked 
by the kink in shock position and an intense burst of the 
LE/ME radiation. From now on, LE and ME luminosities 
begin to oscillate on a time scale roughly corresponding to 
the fundamental mode of the oscillatory instability (Sect. 
3.1). Note that the post-shock region of the oscillating shock 
does not radiate in HE, therefore, no rapid variations of the 
high energy emission are expected. 



For reasons explained in Sect. 3.1 in the next two sec- 
tions we will concentrate our discussion on the spectrum 
emitted by the shock of original velocity v B o = 6000 km s~ x 
(roughly equal to the velocity of the forward shock consid- 
ered by Terlevich et al. 1992). 



4.2 The continuum 

The emergent shock spectrum just after the thin shell for- 
mation (at t ~ 3.20 yr, i.e. at the moment the shock achieves 
its maximum luminosity) is shown in Fig. pi The energy flux 
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Figure 5. Temporal variations of total, HE, ME, and LE luminosity of the post-shock region (thick, medium, thin and dotted lines, 
respectively), (a) initial shock velocity d s o = 2000 km s — 1 ; luminosity scale: 3.71 X 10 7 erg cm - 2 s _1 . (b) initial shock velocity v s o = 6000 
km s _1 ; luminosity scale: 1.52 X 10 s erg cm - 2 s _1 . 



Figure 6. The emergent spectrum of the shock 6000 km s -1 at time t = 3.20 yr, just after the thin shell formation, (a) emergent flux; 
(b) transmitted flux; (c) solid line: sum of the reflected flux and half of the flux generated between the main shock and the shell; dashed: 
reflected flux only. The energy flux vF v (erg cm" 2 s _1 ) is plotted in logarithmic scale. Bottom and top scales arc given in logarithm of 
energy in Rydbergs and keV, respectively. See text for details. 



F v (Fig. ga) is the sum of three components (Figs. ^|(b) and 



(1) the transmitted flux (originating between the shell and the 
left edge of the grid, and partly absorbed by the shell); 

(2) the half of the flux generated between the main shock and 
the shell (emitted directly towards the observer); 

(3) the reflected part of the second half of flux generated be- 
tween the main shock and the shell (originally emitted to- 
wards the shell and reflected by the shell towards the ob- 
server) . 

The emergent spectrum has a power-law shape F v oc v~ a 
with a mean spectral index a m 0.63, and with a high-energy 
cutoff around 10 keV. For energies above few keV the spec- 
trum is completely dominated by the emission produced by 
the hot region between the shell and the left boundary of 
the grid that has not yet had enough time to cool d own due 
to its very high initial temperature (see Sections 3.2 and 



4.1). This part of the spectrum is unlikely to be affected by 
any improvement of cooling function for nonequilibrium ef- 
fects. For lower energies, most of the radiation comes from 
the post-main-shock region located between the main shock 
and the shell. The UV/optical emission is produced by the 
shell ionized predominantly by photons emitted from the 
hot region obscured by the shell. As the emission from this 
region is modeled correctly, the low-energy part of the shock 
spectrum is also calculated correctly. 



4.3 The line emission 

The ionizing flux, as defined in Sect. ^, consists of two com- 
ponents only: 

(1) the flux originating between the shell and the left edge of 
the grid; 

(2) half of the flux generated between the main shock and the 
shell. 
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Table 2. Time lags for selected emission lines. 



line 




time lag [d] 






observed 


TTRFM 


this paper 


H/3 A4861 


20 


30-40 


17.9 


Lya A1216 


12 


1-15 


5.1 


He II A1640 


4-10 


12-20 


5.1 


Mg II A2798 


34-72 


40-50 


40.5 


C III] A1909 


26-32 


20-30 


25.2 


Si IV A1397 


12-34 


12-17 


5.8 


C IV A1549 


8-16 


10-15 


3.3 


N V A1240 


4 


2-10 


5.1 



The ionized medium was identified with the dense shell. The 
shell of the 6000 km s _1 shock contained ~ 50 zones and its 
average thickness was ~4x 10 13 cm. At selected evolution- 
ary times the ionizing spectra together with corresponding 
shell models were used as the input to CLOUDY photoion- 
ization code. 

The output from CLOUDY consisted of intensities and 
intensity ratios of selected lines as functions of time. It was 
found that due to finite time of the shell formation line in- 
tensity maxima lagged behind ionizing flux maximum, the 
lag being shorter for high ionization lines, and longer for 
low ionization ones (Table |). Temporal variations of se- 
lected line intensities and ionizing flux intensities during the 
shell formation phase are shown in Fig. ^. Fig. ^| presents 
temporal variations of selected line ratios throughout the 
simulated evolution together with variations of the total hy- 
drogen column density -/V(H) and maximum density of the 
shell n max . 



5 DISCUSSION 

The general conclusion of our calculations is that the basic 
features of the fast shock model are similar to those observed 
in AGNs. 

First, the emergent spectrum of the shock has a power- 
law shape F v oc v~ a with a mean spectral index of a m 0.63, 
and with a high-energy cutoff around 10 keV. Note that the 
hard X-ray (2 — 10 keV) AGN spectrum is characterized by 
a single power-law with a — 0.7 (Turner & Pounds 1989). 
The high energy end of the spectrum (hereafter the cutoff 
energy) is produced by the gas located behind the shell, and 
we recall that this gas was processed at the beginning of 
evolution by the shock moving with a velocity nearly equal 
to the original shock velocity. Assuming that the cooling was 
ineffective and the temperature of the gas did not change 
appreciately till the time of the shell formation we can state 
that the cutoff energy is directly related to the original shock 
velocity. Therefore, one may expect that the spectrum of the 
shock with the initial velocity equal to the initial velocity of 
a supernova blast wave (v B o = 1 — 2 x 10 4 km s _1 ) would 
have the energy cutoff at around 100 keV, as the post-shock 
temperature is proportional to square of the shock velocity 
(such models were not obtained as we considered evolution 
for t > t sg , see Sect. ^|). This is the cutoff value expected 
for Seyfert galaxies (see Zdziarski, Zycki & Krolik (1993b) 
for references). Note also that our model shock spectrum 
very closely resembles thermal Comptonization component 



used in calculations of X-ray spectra of AGNs (Zdziarski 
et al. 1993b, Zdziarski, Lightman & Maciolek-Niedzwiecki 
1993a). 

Second, intensity ratios of spectral lines produced in the 
photoionized shell agree fairly well with those measured in 
AGNs (Fig. |). The values for AGNs in Fig . | were taken 
after compilation of Kwan & Krolik (1981), and they are 



marked on the vertical axes by dots. We also c onfirm pre- 
liminary results obtained by Terlevich et al. (1992). The 
predicted total hydrogen column densities of the ionized 
medium are of the order of 10 23 cm -2 and the maximum 
gas density is roughly equal to 10 11 cm -3 . Our model fails 
in explanation of low Lya/H/3 ratio, but this problem is 
common for other photoionization models of AGNs. 

Finally, our model exhibits time lag, i.e. a delay between 
maxima of ionizing continuum intensity and line intensities. 
TTRFM using a very simple model of the shell and ioniz- 
ing continuum obtained time lags listed in Table |^. With 
our study, which is entirely free of arbitrary assumptions 
concerning the shell and the ionizing continuum, we cor- 
roborate their findings. The predicted time lags agree well 
with those observed in AGNs. However, this result should be 
treated with some caution as the observational definition of 
the time lag (involving statistical measurements) is different 
from the one employed here. 

Taken together, the above properties provide an indi- 
cation that cSNRs play a significant role in the AGN phe- 
nomenon, as it was proposed in the starburst model (see, 
e.g., Terlevich et al. 1992). One may be almost sure that 
more sophisticated model spectra based on appropriate as- 
semblage of shocks with different velocities would recover 
the ME excess observed in AGNs (i.e. the so called big blue 
bump, Turner & Pounds 1989). This is because in a massive 
starburst there will be many cSNRs coexisting, most of them 
in relatively long lasting advanced evolutionary phases char- 
acterized by low shock velocities. Such models would also 
enable statistical measurements of time lags better suited 
for comparison with observations than these presented here. 

Our conclusion is that cSNR shocks are an attractive 
research field whose exploration has only begun here. The 
finding that time lag phenomena may be related to varia- 
tions of physical parameters of the photoionized medium (as 
originally proposed by TTRFM) seems to be particularly at- 
tractive and it clearly deserves future investigation. 
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Figure 7. Temporal variations of selected line intensities and ionizing flux during formation of the shell for initial shock velocity 
i>sO = 6000 km s _1 . All values are shown in logarithmic scale; time is given in years on horizontal axis. Vertical line marks the time of 
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Figure 8. Temporal variations of selected line ratios and parameters derived from photoionization models for initial shock velocity 
v a o = 6000 km s _1 . Left scale: full line; right scale: dotted line. Line ratios typical for AGNs are marked by dots on vertical axes. All 
values are shown in logarithmic scale; time is given in years on horizontal axis. 
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